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Abstract 

We introduce a class of local stochastic volatility models. Within our framework, we obtain an 
expression for both (i) the price of any European option and (ii) the induced implied volatility smile. To 
illustrate our method, we perform specific computations for a CEV-like model. 
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1 Introduction 



A local volatility model is a stochastic volatility model in which the volatility er t of an asset A is a function 
of the present level of A. That is, a t — a(X t ). Among local stoch astic v olatil ity models, perhaps the most 



well-known is the constant elasticity of variance (CEV) model of ICoxl (|1975f ). An ext ension of the CEV 



model to defaultable assets (the Jump-to-Default CEV or JDCEV model) is derived in 



Carr and Linetskv 



(|2006l ) . One advantage of these two local stochastic volatility models is that they allow for closed-form 
pricing formulas for European options written as infinite scries of special functions. 

In this paper, we introduce a class of local stochastic volatility models which, like the CEV and JDCEV 
models, allow for European option prices to be written down in closed form as an infinite series. Additionally, 
we derive an expression for the exact implied volatility surface induced by our class of models. Previous stud- 
ies of the implied volatility surface induced by local volatility models focused on heat-kernal expansions to de- 



Gatheral. Hsu. Laurence. Ouyang and Wang 



Dupire 



1994) solves the 



rive a s ymptotic approximations of the volatility smile (see e.g 

([20101 ) ; iHenrv-Laborderd (|2005f ) and references therein) . It is worth mentioning that 
inverse problem of finding a formula for the local volatility function the produces a given observed implied 
volatility surface exactly. 
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The rest of this paper proceeds as follows: in section [5] we present our model and assumptions. In section 
[3] we derive a formula for the price of a European option in our modeling framework. In section 3] we provide 
an formula for the implied volatility smile induced by our model. As an example of our framework, in 
section [5] we perform explicit pricing and implied volatility computations for a CEV-like model. Numerical 
results are provided at the conclusion of the text. An appendix with some mathematical background is also 
provided. Concluding remarks can be found in section [5] 



2 Model and Assumptions 

We assume a frictionless market, no arbitrage and take an equivalent martingale measure P chosen by the 
market on a complete filtered probability space (Cl,3 r ,{3 r t,t > 0},P). The filtration {£Ft f i > 0} represents 
the history of the market. All processes defined below live on this space. For simplicity we assume zero 
interest rates and no dividends so that all assets are martingales. We consider a diffusion X with lifetime £ 
whose dynamics are given by 

dX t = (a 2 +e V (\ogX t )) 1 ' 2 X t dWt, (1) 

were, a > 0, e > 0, the function r\ : R — > R + is C§°(R) and W is a Brownian motion Q. Note that X 
has local stochastic volatility <j{Xt) = (a 2 + e 7?(logA t )) 1/2 . Obviously, if 77 = then X is a geometric 
Brownian motion. This will be key for our implied volatility analysis in section 21 Observe that both zero 
and infinity are natural bou n daries according to Feller's boundary classification for one-dimensional diffusions 



(see iBorodin and Salminenl (|2002r ) pp. 14-15). That is, both zero and infinity are unattainable. 

In what follows it will be convenient to introduce Y — log A. A simple application of Ito's formula shows 
that Y satisfies 

dY * = -\ (a 2 +ev(Xt)) dt + ( fl2 + e 17(1*)) 1/2 dW t . 

3 Option Pricing 

We wish to find the time-zero value u e (t,y) of a European-style option with payoff h(Y t ) at time t > 0. 
Using risk-neutral pricing we have 

u e {t,y)=E y h{Y t ), 



1 The notation Cq°(R) indicates the space of infinitely differentiable functions with compact support. 
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where the notation E y indicates expectation starting from y — \ogX$. The function u E (t,y) satisfies the 
Kolmogorov backward equation 

(-d t +A £ )u s = 0, u e (0,y) = h(y). (2) 

where A £ is the generator of the process Y. The domain of A e is defined as the set of / for which the limit 
limt^o j (E y f(Yt) — f(yj) exists in the strong sense. For any / € Cq(R) the generator A e has the explicit 
representation 

A £ =A Q +er l A 1 , A = ^a 2 (d 2 - d) , A 1 = ^ (d 2 - d) , dom(Ai) = Cg(R). 

Remark 1. The operators Aq and Ai are normal operators in the Hilbert space !K = L 2 (R, dy) and satisfy 
the following (improper) eigenvalue equations (neither Aq nor A\ have any proper eigenvalues) 

A fa = Mx, tPx = "4=e iA «, <j> x = \a 2 (-A 2 - iX) , 

V 27T ^ 

Ai1>\ = x\1>\, fa = ^=e x \ xx = \ (-a 2 - iX) . 

V Z7T ^ 

Note that the eigenfunctions satisfy ip^) = f tpx[y)ipfi(y)dy — 8{X — /i). Note also that Borel-measurable 
functions of normal operators (e.g., g(Ao)) are well-defined, as explained in Appendix 1X1 

We seek a solution to Cauchy problem ([2]) of the form 

oo 

= 5>™u„. (3) 



OO 
71=0 



We will justify this expansion in Theorem [2] Inserting the expansion ([3]) into Cauchy problem ([2]) and 
collecting terms of like powers of e we obtain 

0(1) : (-d t + A )u = 0, u (0, y) = h(y), 

0(e n ) : (-d t +A )u n = -rjAiUn-i, u n (0,y)=0. 

The solution to the above equations is 

0(1): u (t,y) = e tA °h(y), 

0(e n ): u n (t,y)= [ dse^-^^AiUn-x^y). 

Jo 

Using the equation from appendix [Al we obtain 

0(1): u (t,y)= [ dXe^^h^xiy), 

JR 

0(e"): u n (t,y)= [ [ dsdfie^** (^,^i^n-i(s, •)) ^(v), 
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After a bit of algebra, we find an explicit representation for u n (t,y) 

n+l 

We have now obtained a formal expansion for the price of a European option. The following theorem provides 
conditions under which the expansion is guaranteed to be valid. 



Theorem 2 (Option Price). Suppose e < t^it, where \\ri\\ = \J (77, rj). Then the option price u e (t,y) is given 
by © - ©. 

Proof. See Appendix IB"1 □ 



4 Implied Volatility 

In this section we fix (t, y) and a call option payoff h(y) = (e y — e k ) + . Note that 

p k—ik\ 

(^»--=— -— , Im(A)<-l. 
The following definitions will be useful: 

Definition 3. The Black-Scholes Price u BS : R+ -S- R+ is defined as 



(a) := f dXe^W^h)^, ^ s (a) = l -a 2 {-\ 2 - iX). 



Definition 4. The Implied Volatility is defined implicitly as the unique number o~ e € M + such that 

u BS {a e ) = u £ , (5) 

where u £ is as given in Theorem [2J 



Remark 5. Note that uq = u BS (a). As shown in iLoria (|2012l ). when u £ can be expanded as a power series 
whose first term corresponds to u BS , one can obtain the exact implied volatility surface corresponding to 
u e . 

Remark 6. For < t < 00 the existence and uniqueness of the implied volatility a e can be deduced by 
using the general arbitrage bounds for call prices and the monotonicity of u B . 

Remark 7. Note that u BS is an invertible analytic function that satisfies d p u BS (a) > for all a > 0. By 
the Lagrange inversion theorem, the inverse [u B ] _1 of such a function is also analytic. 
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Clearly, u e is an analytic function of e (we derived its po wer series expansion). It is a useful fact that the 
composition of two analytic functions is also analytic (see iBrown and Churchilll (|1996T ). section 24, p. 74). 
Thus, in light of Remark [71 we deduce that a e = [u ss ] _1 (u' r ) is an analytic function and therefore has a 
power series expansion in s. We write this expansion as follows 



<r e = a + <5 £ , 

Taylor expanding u BS about the point do we have 
u BS (^)=u BS (a»+n 

oo 1 



CO 



(6) 



fe=i 



n=0 

u BS 

„,BS 



oo / oo \ 71 

n=l ' \k=l / 



u Bb (a Q ) 



oo 1 

n=l 
oo 

M + E e * 



e e n 

fc=l VilH hi„=fc i=l 



fc=i 



9> BS (a ) 



,BS( 



oo 1 / n \ 

"=1 V'i+--+3„=fei=l / 

^+e^( e n**W 

n=2 yi+-+J»=fci=l / 



; K) + E £ 

fc=l 

Now, we insert expansions © and (UJ into ([5]) and collect terms of like order in e 
0(1) : u = u BS {a a ), 

0(e k ) : u k 

Solving the above equations for {ofc}^L we find 
0(1): a = a, 

0{e k ) : 



u Bb (a ). 



(7) 



n=2 " \ji+-+j„=fci=l / 



k > 1. 



Er^f E lUW M («*>) 

n=2 ' \ji+-+j„=fci=l / / 



fc > 1. (8) 



Remark 8. The right hand side of ([5]) involves only cr,- for j < fc — 1. Thus, the {cfc}^^ can be found 
recursively. 

Remark 9. Note that d™u BS ((j) is easily computed using 
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Explicitly, up to 0(e 4 ) we have 



d a u a ' 

1 2 n2 

d a u 

u 3 - (<j 2 cndl + ji<jfdl)u 
d a u 

u A - {ozoxdl + \o\dl + \o-io\dl + ±<rfd*)u 
d a u 

We summarize our implied volatility result in the following theorem: 

Theorem 10 (Implied Volatility). The implied volatility a e defined in (O is given explicitly by (J6j) where 
Co = a and {crk}^ =1 are given by (JSj) . 

Remark 11. Everything we have done so far is exact . The accuracy of the implied volatility expansion ([6]) 
is limited only by the number of terms one wishes to compute. 



0(e) : a x = 

0(e 2 ) : o- 2 = 

0(e 3 ) : a 3 = 



Coxl (|l975[ ) the dynamics of X are assumed to be of the 



5 CEV-like Example 

In the constant elasticity of variance (CEV) model of 
form dXt = •JeX^' 2 X t dWt- A key feature of the CEV model is that, when j3 < 0, volatility a{x) = ^fex 13 / 2 
increases as x \ 0, which (i) is consistent with the leverage effect and (ii) results in a negative implied 
volatility skew. However, values of j3 < also cause the volatility to drop unrealistically close to zero as x 
increases. If we choose rj(y) = e@(y) := e@ v in (Q]) then the dynamics of X become 

dX t = (a 2 +eXP) 1 ' 2 X t dW u 

Note that the local volatility function o~(x) — (a 2 + ex^) 1 / 2 behaves like a(x) ~ ^fex^l 2 as x \ and 
behaves like a constant o~(x) ~oasi /* oo. See figure [1] for a comparison of our local volatility function 
and the CEV local volatility function. 

Remark 12. Because e^ y is unbounded as y —> — oo (recall (3 < 0), the function ep ^ Cp°(R). However, 
we can modify the domain of u £ to be R + x M where Mo := (j/0i °°)- The operators Aq an d -A-i would then 
be defined on L 2 (M. Q ,dy) and the domain of these operators would include a absorbing boundary condition 
at yo (signifying default of X the first time X reaches the level e Vo ). Note that || e^a || := (f°^ lep^dy) 1 / 2 = 
e^ va J \J — 2(3. In the analysis that follows, it will simplifyy computations considerably if we continue to work on 
L 2 (R, dy) as working on L 2 (Rq, dy) would require modifying the eigenfunctions ip\ from complex exponentials 



G 



exp(iAy) to cosines cos(Xy). However, the simplification comes at a cost; in light of the conditions of theorem 
© our results may not be valid for values of y < — ^ log a ^~ 2/3 ■ 

We wish to find a simplified expression for u n (j4|) for the case r) = ep. Noting that 

epAitpx) = Xa 5(\ — (i— i(3), 



dX \Z^ TT» (A 7. \ 11 Oa, h) ipx-inp 



we see that the n + 1-fold integral (|4| collapses into a single integral 

= e n p / d\ \ II XA-ifc/3 (ipx,h)ijjx- (9) 

Remark 13. Although we have written the option price as an infinite series ([3]), from a practical standpoint , 
one may only compute u e ~ S^Lo £Tl-u n- F° r an y finite iV we may pass the sum through the integral 
appearing in @. Thus, for the purposes of computation, the best way express the approximate option price 
is 

u e k d\ (ipx,h) VA V £ n e n/3 V TT xx-ikp ■ 

Note, to obtain the approximate value of u £ , onZy a single integration is required. This makes our pricing 
formula as efficient as other models in which option prices are expressed as a Fourier-type integral (e.g. Levy 
processes, Heston model, etc.). 

Some Numerical Results 

Define the transition density p s (t, y; y ) and the 0(e n ) approximation of the transition density p^[t, y; yg) 

n 

P e (t,y;y Q )= E yo 5 y {Y t ) , p^ (t, y; y ) = ^ ^Vfc (*, V, tfo) • 

In figure^] we plot the approximate transition density p( n ). Next, define the 0(e n ) approximation of the 
implied volatility 

n 
k=0 

where the Ok are given by ([5]) . In figure [3] we provide a numerical example illustrating convergence of 
er(™) to a e . We compute a e by calculating u e first using Theorem [5] and then by inverting the Black- 
Scholes formula numerically. We plot implied volatility as a function of the log-moneyness to maturity ratio, 
LMMR := (k — y)/t. Convergence is fastest for values of k near y and slows as k moves away from y. 
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6 Conclusion 



In this paper we introduce a class of local stochastic volatility models. Within our modeling framework, 
we obtain a formula (written as an infinite series) for the price of any European option. Additionally, we 
obtain an explicit expression for the implied volatility smile induced by our class of models. As an example 
of our framework, we introduce a CEV-like model, which corrects one possible short-coming of the CEV 
model; namely, our choice of volatility does not drop to zero as the value of the underlying increases. In the 
CEV-like framework, we show that option prices can be computed with the same level of efficiency as other 
models in which option prices are computed as Fourier-type integrals. 

Thanks 

The author would like to thank Bjorn Birnir for his helpful comments. 
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A Spectral theory of normal operators in a Hilbert space 



In this appendix we summarize the theory of normal o perators acting; on a H ilbert 



Reed and Simonl (|1980f ) and 



space. A d etailed 



Rudinl ( 1973 ) 



exposition on this topic (including proofs) can be found in 

Let !K be a Hilbert space with inner product (•, •). A linear operator is a pair (dom(A)^A) where dom(A) 
is a linear subset of J£ and A is a linear map A : dom(A) — > 5£. The adjoint of an operator A is an operator 
A* such that (Af,g) = (f,A*g),Vfe dom(A),g € dom(.A*), where 

dom(A*) :={ 9 €J(:3fteJ{ such that (Af,g) = (f,h)Vf € dom(yi)}. 

An operator (dom(A),A) is said to be self-adjoint in 3-C if 

dom(A) = dom(A* ) , (Af, g) = (f,Ag) V /, g e dom(A) . 

Throughout this appendix, for any self-adjoint opera tor A, we will assume that dom(A) is a dense subset of r K. 



A densely defined self-adjoint operator is closed (see 



Rudinl (|1973f ). Theorem 13.9). An operator (dom(A),.A) 



is said to be normal in Jt if it is closed, densely defined and commutes with its adjoint: A* A = AA* . Clearly, 
every self-adjoint operator is a normal operator. 

Given a linear operator A, the resolvent set p(A) is defined as the set of A € C such that the mapping 
(A — Id A) is one-to-one and R\ := (A — Id A) -1 is continuous with dom(i?>) = JC. The operator R\ : — >• 5£ 
is called the resolvent. The spectrum a (A) of an operator A is defined as a(A) := C \ p{A). We say that 
A € cr(.A) is an eigenvalue A if there exists iji G dom(A) such that the eigenvalue equation is satisfied 



(10) 



A function -0 that solves (|10p is called an eigenfunction of .A corresponding to A. The multiplicity of an 
eigenvalue A is the number of linearly independent eigenfunctions for which equation (|10[) is satisfied. The 
spectrum of an operator A can be decomposed into two disjoint sets called the discrete and essential [^| 
spectra: a(A) = ad(A) U a e (A). For a normal operator A, a number A € C belongs too^L/l) if and only if A 



Rudinl (|l973l) . Theorem 12.29). 



is an isolated point of a(A) and A is an eigenvalue of finite multiplicity (see 

A projection-valued measure on the measure space (C, 23(C)) is a family of bounded linear operators 
{E(B),B £ 3(C)} in J£ that satisfies: 

1. E(0) = and E(C) = Id. 

2. is an orthogonal projection. That is, E 2 (B) = E(B) and E(B) is self-adjoint: E*(B) = E(B). 



2 The essential spectrum may be further decomposed into the continuous spec trum and the residual spectrum. It can be 
shown that the residual spectrum of an ordinary differential operator is empty (see 



Roach 



(1982), page 184). 
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3. E(Af]B) = E(A)E{B). 

4. If B = B t and B l n B 3 ■ = for i ^ j then = lim^oo £ ™=i ^(^0. where the limit is in thc 
strong operator topology. 

5. For every in /, g £ "K the set function fj,t !g (B) := (f,E(B)g) is a complex measure on 23(C). 

Theorem 14 (Spectral Representation Theorem). There is a one-to-one correspondence between normal 
operators A and projection-valued measures E on "K, the correspondence being given by 



A= / XE(dX). 



If * s a B orel function on C then 



g(A) = I g(X)E(dX), 



Proof. See 



domfoGA)) = {/ € JC : / |ff(A)|V/,/(^) < oo}. (11) 

«/<r(.A) 



Rudinl (|l973l ) Theorems 12.21 and 13.33. 



□ 



As a practical matter, if .A is an differential operator acting on a Hilbcrt space L 2 (M.,dy), then the operators 
defined by (|ll|l can be constructed by solving thc proper and improper eigenvalue problems 



proper: 
improper: 



A ip n = (j) n ijj n , 
A ip\ — cj)\ tpx, 



4>n € (Td(A), 

4>x e a e {A) 1 



For the improper eigenvalue problem one extends the domain of A to include functions all functions tp for 
which Af makes sense and for which the following boundedness conditions are satisfied 

lim \ip{y)\ 2 < oo. 

y— >±oo 

After normalizing, the proper and improper eigenfunctions A satisfy the following orthogonality relations 

(lpn,i>m) = S n , m , (lp\,lpy) = 8(X - A'), (lp n ,1p\) = 0. 



The operator g(A) in (fTTj) is constructed as follows (see lHanson and Yakovlevl (|2002h . section 5.3.2) 



g(A)f= V g(<f>x)(ipx,f)i>x+ [ g{<t>x){^xJ)i>\dX 



\£a d (A) 



It is not always easy to evaluate divergent integrals of the form (ipXi^X 1 ) an d verify that they are in fact 
delta functions S(X — A') . A method for d irectly obtaining properly normalised improper eigenfunctions can 
be found on page 238 of iFriedmanl ()1956l ). 



3 The term "improper" is used because the improper eigenvalues A ^ od{A) and the improper eigenfunctions ipx t since 
(i>X,">P\) = co- 
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B Proof of Theorem [2] 

Our strategy is to show that A £ = Aq + er\A\ generates a semigroup 7f = exp(tA £ ). This will guarantee 
that u £ = y e t h(y) is an analytic function of e, which in turn, justifies the use of expansion Throughout 
this section we will work on the Hilbert space J£ = L 2 (R, dy). We let do m(./lj) = Cn°(R) and we note that 



Chernoffl ( 1972 ) 



(R) is a dense subset of 3-C. Our analysis begins with a Theorem from 

Theorem 15. Let A be the generator of a Co contraction semigroup = cxp(tA) on a Banach space. Let 
£ 25 be a dissipative operator with a densely defined adjoint. Assume that the inequality 

WeftuW < a \\u\\ +b\\Au\\ , VtiG dom(A) 

holds for some a > and b < 1 (i.e., the operator sH is A-bounded with bound b < 1). Then the closure of 
A £ := A + eH generates a Cq contraction semigroup T £ = exp(L4 e ). 

Remark 16. The operator ^lo is the generator of a Co contraction semigroup CP^ = exp(tAo) on !K. 
To show that er\A\ is dissipative, the following Theorem will be useful: 

Theorem 17. Let A be a linear operator with domain dom(A) = C 3O (M). Then A satisfies the positive 
maximum principle if and only if 

1 O , , „o . , , „ f , . , / g 



A = 2« 2 (y)9 2 + Ky)d + J v{y, dz) {e zd - 1 - I {z<R} zd) - c(y), (12) 
for some a{x) > 0, b(x) € R, c(x) >0,R€ [0, oo] and v{y,dz) satisfying 

v{y,dz) (l A z 2 ) < oo. 
Operators of the form (|12j) are called Levy- type operators. 

Proof. See Theorem 2.12 of|HobJ |l998l ). □ 



Rem ark 18. An operator that satisfies the positive maximum principle is dissipative (see 



Ethier and Kurtz 



(|1986l ). Lemma 4.2.1 on page 165). 



Remark 19. The operator er\A\ is clearly of the form (IT21 . Therefore, er\A\ is dissipative. 

Remark 20. The adjoint of er/Ai, given by (eijAi)* = eA*rj, has domain dom(eAlr]) — C X3 (R), and is 
therefore densely defined in "K. 

Since Aq generates a Co semigroup and er\A\ is dissipative and has a densely defined adjoint, we have only 
to show that et]A\ is .Ao-bounded with b < 1. 
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Proposition 21. Suppose r\ <G Cq°(M) and e < t^tt (which is the condition given in Theorem^. Then 
srjAi is Ag-bounded with a = and b < 1. 

Proof. Clearly, for any u 6 dom(.Ao) we have 

H^yix U || < e \\n\\ ■ \\a iU \\ = 4 IMI ■ ll^ll ^ IKo^ll . 

□ 

The proof of Theorem [2] is complete. 
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Figure 1: A comparison of the CEV volatility a{e y ) = ^/ee^ y / 2 (dashed black) and our "shifted" version 
of the CEV volatility a{e y ) = (a 2 + ee^) 1 / 2 (solid blue). Notice that the CEV volatility drops to zero as 
y — > oo whereas our shifted version stays above a. The following parameters are used in this plot: a = 0.20, 
Ve = 0.15, /3 = -0.75. 
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Figure 2: A plot of the approximate transition density p( n > [t, y; 0) (solid blue) for different values of n. For 
comparison, we also plot p^ (dashed black). Note that the density of Y± has a fat tail to the left, which 
is expected since a(e y ) increase as y — > — oo. The following parameters are used in these plots: a = 0.20, 
= 0.15, /3 = -0.75, t = 2.0. 
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